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I. INTRODUCTION 



The Higgs boson is the only missing, undiscovered component of the Standard Model of 
Particle Physics. Within the context of the Standard Model, the Higgs boson is responsible 
for the breaking of the SUl{2) x U{1) gauge symmetry and provides the mechanism for the 
generation of masses of the corresponding gauge bosons: and Z. In addition, the same 
mechanism provides masses for the leptons and quarks via Yukawa couplings. Therefore, 
the discovery and subsequent study of the Higgs boson properties is of the highest priority 
for particles physics community. 

A precise theoretical understanding of Higgs production rate is critical to any attempts 
to search for the particle. The dominant production mechanism for Higgs bosons in hadron- 
hadron colliders is via gluon-gluon fusion pp ^ gg ^ H, in which gluons fuse through 
a virtual top quark triangle to produce the Higgs. Such a process provides the largest 
production rate for the entire Higgs mass range of interest. 

The standard description of hard processes in hadron collisions within the framework 
of QCD parton model, which reduces the hadron-hadron interactions to the parton-parton 
ones via the formalism of the hadron structure functions. The most popular approach is the 
QCD coUinear approximation based on the well known collinear factorization theorem 
[3|. In this approach all particles involved are assumed to be on the mass shell, carrying only 
longitudinal momenta, and the cross section is averaged over two transverse polarizations 
of the incident gluons. The transverse momenta of the incident partons are neglected in 
the QCD matrix elements. However, at small x, the effects of finite transverse momenta of 
partons become increasingly important, especially in the analysis of jets and heavy-quark 
production, and there is no reason to neglect the transverse momenta of the gluons in 
comparison to the quark mass. The method to incorporate the incident parton transverse 
momenta is referred to as /cj_-factorization approach [4, 5]. Here the Feynman diagrams are 
calculated taking into account the virtualities and all possible polarizations of the incident 
partons. There are widely discussed applications of the /c_L-factorization approach to hard 
QCD processes like the J/ip hadroproduction |6^, charmonium production heavy quark 
photo- ,9|] and hadroproduction [lOj, Higgs boson production 11, l3|. Some exclusive 



processes in the framework of the fc^-factorization approach are described in detail in [Ic 

In the lowest order the gluon coupling to the Higgs boson in the Standard Model is 
mediated by triangular loops with top quarks as shown by the Feynman diagram in Fig. ^ 
This process for on-shell gluons p^, fc^ = is well-known in the literature 114], so we will 




FIG. 1: Gluon coupling to the Higgs boson via i-quark triangle loop. 
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For a large quark mass r ^ 1 the form factor Ji can be expanded in powers of r ^ 
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The form factor I2 is not relevant in the coUinear approximation because it does not come 
into the squared matrix element in the on-shell limit. 

In order to use the fc^-factorization approach we will take into account the non-trivial 
virtualities of external gluons. The main technical problem here is to calculate the vector 
and second-rank tensor Feynman integrals with p^^k"^ 7^ 0. But the symmetric properties of 
full amplitude together with a usage of convenient projectors allow to avoid such difficulties 
and represent all relevant form factors in terms of scalar three- and two-point functions plus 
extra finite subtractions appearing due to the regularization. 

In the present work we shall study in detail the fusion of off-shell gluons in proton-proton 
collisions producing a scalar Higgs boson by the top quark triangular loops with finite (not 
infinite) top quark mass. We analytically calculate the exact amplitude of this process in 
terms of two relevant form factors. The amplitude is studied in various kinematical regions. 
Furthermore we discuss different fermion mass limits. We find that the squared matrix 
element taking into account the non-zeroth gluon virtualities slightly (by several percents) 
increases with the growth of gluon transverse momentum. We also discuss consequences of 
the non-zeroth virtualities for scalar boson production in the /c_L-factorization approach for 
large transverse momenta. 



II. THE AMPLITUDE FOR THE HIGGS BOSON PRODUCTION IN A FUSION 
OF OFF-SHELL GLUONS 

Let us start from a general tensor representation of the triangle amplitude shown in Fig.^ 

T^,y{k, p) = g^^Fi + kf,k^F2 + Pf^Pi^F^ + k^^p^F^ + p^k^F^. (2.1) 

Here Fj = Fj{q'^,k'^,p'^,m'j), j = 1,...,5 are the Lorentz invariant form factors, is 
the fermion mass in the loop. The Bose symmetry of the amplitude T^^{k,p) = T^^{p, k) is 
equivalent to 

F,{k,p) = Fi(p,fc), F2{k,p) = Fs{p,k), 

F,ik,p) = F,ip,k), F,ik,p) = F,{p,k). (2.2) 

The gauge invariance leads to the vector Ward identities p^T^j^^ = 0, k'^T^^, = which in 
terms of form factors gives 

Fi + /F3 + {kp)F^ = 0, {kp)F2 + p^F, = 0, 

Fi + eF2 + {kp)F^ = 0, {kp)F3 + k^F^ = 0. (2.3) 
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Therefore we have finally only three relations 

F,+p^F3 + {kp)F4 = , 
{kp)F2+p^F, = 0, 
p^Fs - k^F2 = . 



(2.4) 



This reduces the number of independent form factors to two. So the tensor representation 
of the amplitude can be expressed in the following form 



T^,u{k,p) = i6 
+ 
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where is the strong coupling constant, v = {Gf\/2)^^^'^ is the electro-weak mass scale, a, b 
are the color indices for two off-shell gluons. This general form of the amplitude coincides 
with the standard expression for on-shell gluons with 
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In order to calculate the form factors Gi, G2 we introduce two projectors. The most 
convenient and simple choice is one symmetric P^^^ = g^u and one antisymmetric Q/iv = 
Pfiku — Pvk^ 2-rank tensors. Corresponding projections are 



^1 = T.^Q''' = z6^'^- ( [kY - {kpf]G^ + 
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Solving this system of equation with respect to Gi, G2 we get 

^,,,_,nv {kp)S, + {{kp)^-kY)S2 
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(2.7) 
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The explicit expression for the amplitude reads 
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i{2T[)^ [(r — k)"^ — m'j] [r^ — mj-] [{r + pY" 
M^^{k,p) = Tr[(m/+ /- /;)7^(m+ /)7^(m/+ /+ ^)] . 
Employing now the dimensional regularization and using the Passarino-Veltman reduction 
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we get 



^1 = ^5^'— ^ [Coikp)ie+p' + 2(M) - ik' + {kp))Bl' + 



(2.8) 
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where n = 4 — e is the dimension of the space, B^"^, B^^, B^^ are the scalar two-point 
functions 
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and Co is the convergent scalar three-point function. Following Ref. 



i(27r)'' [(r — ky — m'j] [r^ — m'j] [{r + pY — mj] 
= >c{k\p\q') + x{q\k\p') + K{p\q\k'). 
Above for simplicity we have introduced the following notations 
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Of course, the projections Si and S2 are convergent in the limit e and then 
= z6-'^ml [Co{kp){k^ + p^ + 2{kp)) + {{kp) + k'')Li - 

7TV ■' 



So 



ik^+p^ + 2ikp))L, + iikp)+p^)L,], 



id"^— m] [Co(4m2 - fc^ - _ 4(^p)) - Li + 2L2 - ^3 + 2]. 



(2.10) 
(2.11) 
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There is a finite subtraction term in 5*2 as a consequence of regularization. Substituting 
tfiese expressions into ()2.7|) we finally obtain 



2 

TTLf 

_ I ^ 

X [{Am]{{kpf - kY) - 2{kp){{kpf - 2kY) + fcV(fc' + p'))Co + 

+ eip" + {kp))u - {2ky + {kp){e +p^))L2 + 

+ p\e + {kp))L, + 2{{kpf - fcV)], (2.12) 
{{kpY — k'^p^Y 

X [(4mJ((A;p)2 - fcV) + (A;^ + p'){2{kpf + fcV) + 2{kp){{kpf + 2ky))Co + 

+ {2{kpf + ?,{kp)e + A;V)Li - (3(M(^' + 2(2(M=^ + A;V))L2 + 

+ {2{kpf + 3(A;py + fcV)L3 + 2((A;p)2 - fcV)]. (2.13) 

Now the off-shell amplitude can be calculated as M = T^,^A;^j_fc2_L/|ki_L| |k2_L| with T^^, 
given by Eq. (j2.5|) and Gi and G2 given by Eqs. (j2.12j) and (j2.13|) . respectively. 



III. ESTIMATING THE SIZE OF THE OFF-SHELL EFFECTS 



In the so-called heavy fermion limit ^ 00 we get 
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(3.1) 



for gluons with arbitrary virtualities k"^ and p"^. The first limit for Gi coincides with the 
classical on-shell result (|1.3p . 

It is easy to show also that the form factor Gi in the limit k'^,p'^ —>■ 
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that coincides with Ji fll.2|) for the relevant case of Higgs boson production from on-shell 
gluon fusion with = (A; + p)"^ ~ m^, ruf = nitop- In this limit the second form factor G2 
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Expansions for r ^ 1 and p^, /c^ = have the following form 
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Let's write analogous expansions with taking into account the non-zeroth gluon virtualities. 

The infinitely heavy fermion limits ()3.1|) do not contain at all the dimensional quantities 
such as g^, p^, /c^, so there is no difference in the order of limits: at first — > and then 
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or vice versa. Therefore it is more convenient to work with the dimensionless 
parameters defined as 



X 



Amj 



< 0. 



(3.5) 



Then the on-shell hmit p^, A;^ ^ is equivalent to ^, ?7 0. The heavy quark hmit ()3.ip 
corresponds to Xi V ^ 0- We can now take into account the terms of first order in ^ and 
7] in ()3.4j) . In the case of Higgs boson production on average m\ ^ K^^)|, so we 

have to take into account the powers of x higher than the powers of ^, t]. Expansions of Gi 
and G2 in up to second order in x and first order in ^ and t] give 
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(3.6) 



We have compared the exact form factors ()2.12|) and ()2.13|) with their expansion counterparts 
()3.6|) for reahstic parameters = 0.17 TeV, rrih = 0.15 TeV. In conclusion, we can use 
the form factor expansions up to |.^|, |?7| = 0.3 with a maximal error of less than 1%. 

To make our calculations useful for the case of Higgs boson production let us turn to the 
physical parameters relevant for proton-proton collisions. Introducing p = Xipi + kij_, k = 
X2P2+k2j_, where ki±, k2± are space-like four- vectors associated with the transverse momenta 
of gluons kij_ and ki_L, pi, p2 are the hadron momenta, so {P1P2) = s/2, and neglecting 



the hadron mass nip -C rrih, rritop, we see that p 



-k2 



< and k"^ 



-k2 



< 0. 



The subject of our analysis is the following normalized projection of the amplitude M 
k2±|, given by the formula 



Tij,i/k^j_k2j_/\^i± 
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COS( 



Gi - (3.7) 
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where Gi, G2 can be taken from ()3.6p for not too large |ki_|_| and |k2_L|, is the azimuthal 
angle between gluon transverse momenta ki_|_ and k2_L, the transverse momentum of the 
produced Higgs boson is = ki^ + k2_L and the virtual gluon polarization tensor has been 
taken in the form ^| 



Aj Ay J 

"ki" 



For brevity, we shall denote the first term of the amplitude by A^i and the second term by 
M.2- Neglecting the second term in ()3.7|1 and taking the on-shell value of Gi one obtains 
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(3.8) 



Gl 



This coincides with the amplitude obtained in Ref. |11| within a numerical factor. 
Let us now quantify some of the off-shell effects discussed above. 
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A. Effect on form factors 



Let us start with a simple case of form factors. In Fig. |21we show the dependence of the 
two off-shell form factors Gi and G2 on parameters ^ and rj. The results are normalized to 
the on-shell values of the form factors 
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In this calculation mtop=0.17 TeV, 777,^=0.15 TeV. We see that the first form factor Gi 
slightly drops, while the second form factor G2 grows with increasing |^| and \ri\. The first 
form factor with taking into account of finite quark mass differs from one obtained in the 
framework of effective approach by 5 ' 



70. 



B. Effect on amplitude 



By averaging the amplitude squared over (f) we obtain 
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One can see that in the collinear limit ki^ 



+ k2 ^2 



G1G2 



^2±J 
k2± 



the averaged square of matrix ele- 
ment (M^)(^ coincides with squared matrix element in the covariant gauge T^^T^ ^ g^fi'duv' 
multiplied by 2. 

Let's define the matrix element Mq that is obtained from M given by ()3.7|) by substituting 
the off-shell form factors by the on-shell ones G\ and G^ (see Eq. ()3.10|) ). 



Mn 
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(3.11) 



In Fig. ISlwe show (M^)^/ (Mq)^ as a function of the reduced parameters ^ and rj. One can 
see a drop of (M^)^ relative to (Mq)^ due to the drop of the form factor Gi. 

Let us first estimate the off-shell effect. Following Ref. we take the values of transverse 
gluon momenta in the interval |ki^2±| ~ 5 — 50 GeV. Then the relative effect of replacing 
Gl and G2 by their on-shell counterparts G? and G2 is 



M)<^ - (M^), {G\f - Gl 



(G?)^ 



0.0003 - 0.03 



(3.12) 



correspondingly for the lower and upper limits of gluon transverse momenta. A rather small 
drop of the averaged squared matrix element can be observed, mainly due to the drop of 
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the first form factor Gi. The effect is of the order of 1% or less at typical gluon transverse 
momenta. Therefore the effect of non-zeroth gluon virtualities in the form factors on the 
averaged squared matrix element is rather small. 

It seems interesting to study the behavior of the matrix element in this part of the 
phase space where the first term is small. Actually, when cos0 = the contribution of the 
form factor Gi disappears from (j3.7p . and we have the following simple expression for the 
amplitude 



M\^=^/2 = ^5'^'^-|ki±l|k2±| G2, (3.13) 
In V 



that is determined exclusively by the second form factor 6*2 • As a result there are significant 
consequences of the non-zeroth virtualities on the angular distribution around = 7r/2. 
Then 1,^=^/2 grows with increasing |^| and \ri\. The relative growth is determined com- 
pletely by the growth of the form factor G2 

2 _ qo"^ 

' ' = 0.004- 0.44, (3.14) 
^2 

at |ki 2_L| ~ 5 — 50 GeV respectively. Thus taking into account the non-zeroth virtualities 
could considerably change the angular distribution at (p ^ tt /2. If this can be observed in 
reality it will be discussed in the next section. 
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C. Effect on inclusive cross section 

Let us try to evaluate how big can be the observable effect. This requires a convolution 
of the off-shell subprocess cross section with realistic unintegrated gluon distributions. The 
inclusive cross section for the production of the Higgs boson in hadron-hadron collisions can 
be written as 



dxi dx2 (i^ki_L (i^k2_L 
a = I 

Xi X2 IT IT 
(To//-she«Kkl±,k2±) A{Xi,klj_, flp) A{X2, klj_, flp) , 



(3.15) 



where A{x, k^,/i^) are the unintegrated gluon distributions in the colliding hadrons hi and 

h2. 

Now the distribution in azimuthal angle (p between ki_|_ and k2_L can be calculated as: 



dxi dx2 1 

Xi X2 TT^ 
2 ,,2\ 



ki±dkn^k2±dk2 



^of f— shell (kl±,k2_L) ^(Xi,k^^,/i ) ^(x2,k2_L,/i 



(3.16) 



Inserting the off-shell matrix element squared to the formula for the off-shell subprocess 
cross section we get 
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256 TT^ v'^xiX2sml 



cos Gi 
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(3.17) 



X A{xi,klj^, (i'^)A{x2,kl^, fi'^) dkjj^dkl^dyH, 
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where yn is the center-of-mass Higgs boson rapidity and the longitudinal momentum frac- 
tions must be evaluated as 



3^1 = y ^ ^ ^ exp(yH), X2 = y ^ ^ ^ exp(-yH). 

Expression ()3.17j) . in the on-shell limit for form factors, is consistent with the analogous 
formula obtained by Lipatov and Zotov in Ref. 11 . 

The second term of the off-shell amphtude ()3.7p is new and was not discussed so far in 
the literature. For illustrating the role of the second part of the amplitude in Fig. |3 we 
show the azimuthal angle distributions. In this calculation the BFKL unintegrated gluon 
distributions were used for example (for more details see e.g. There is only a small 

difference between the result obtained with the sum of both amplitudes (A^i + A^2) and 
the result obtained with the first term { M.i) only. The difference becomes visible only 
close to = 7r/2, i.e. when the first amplitude vanishes. We show in the figure also the 
modulus of the interference term. The latter is much smaller than the cross section of the 
first term only, except of ~ 7r/2. Because this effect is negligible except of extremely close 
to = 7r/2, we see no easy way to identify the off-shell effects. 

From the theoretical point of view the exact off-shell matrix element would break the 
familiar cos^ 0-dependence of the cross section. The second term of our exact matrix element 
leads to an asymmetry around 7r/2: 

do /'■K \ da fn \ ,^ 

-<^o <37 o+<^o (3.18) 



2 d(t)\2 

for 00 > 0. In Fig. El we show the ratio of the interference term to the sum of the two 
noninterference terms. Here the off-shell effect is only at the 10~^ level. A bigger effect is 
obtained at ~ 7r/2. A change of sign at = 7r/2 is clearly visible. Putting an extra cut 
on |p_l| could increase the relative asymmetry. In Fig. El we show the ratio obtained with 
such an extra cut |p_|_| > 50 GeV by the dash-dotted line. 



IV. CONCLUSION 

In the present work we have analyzed the effect of the non-zeroth virtualities of external 
gluons on the amplitude of a scalar Higgs boson production. An off-shell matrix element 
for this process was calculated. We have found a new term in the amplitude compared to a 
recent effective Lagrangian calculation. 

A straightforward application of our analysis is in the case of inclusive Higgs boson 
production in proton-proton collisions at LHC We have estimated that the relative drop 
of the averaged square of the matrix element caused by replacing the on-shell form factors 
by the off-shell ones is only about 1% or less at relevant physical parameters of the future 
experiments, so this effect could be verified in the high-precision experiments only. However, 
the effect of taking into account the non-zeroth virtualities on the angular distribution at 
^ 7r/2 is much more significant due to a quick growth of the second form factor G2 as a 
function of gluon transverse momenta. The relative growth of the squared matrix element 
at = 7r/2 is up to 45 % for |ki,2±| ~ 50 GeV. 

The observable effect was estimated in the framework of the /c^-factorization approach by 
convoluting the squared off-shell matrix element and the unintegrated gluon distributions. 
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The effect found is, however, extremely smah and concentrated around (f) — 7r/2. This will 
be extremely difficult to identify in the future rather low-statistics experiments. 

We have discussed deviations from the cos^ 0-dependence of the cross section. A small 
asymmetry around 4> = 7r/2, at the 10~^ level, was found. Taking into account that the 
0-dependence is not directly a measurable quantity the predicted effect seems extremely 
difficult to be identified experimentally. 



We are grateful to A.E. Dorokhov, I.V. Anikin and J. Horejsi for useful discussions and 
comments. This work was partially supported by the grant RFBR 03-02-16816 and by the 
grant of the Polish Ministry of Scientific Research and Information Technology number 1 
P03B 028 28. 
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FIG. 2: Off-shell form factors Gi (left) and G2 (right) normalized to their on-shell values, as 



functions of ^ = — 




FIG. 3: Averaged square of off-shell matrix element {M'^)^ normalized to its on-shell value (Mq)^ 
as a function of = — k^_|_/4mj and t] = — k|^/4m^ . 
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FIG. 4: Azimuthal angle distribution of the cross section. In this calculation the BFKL UGDF 
was used and —2 < yu < 2. The solid line represents the calculation with the full amplitude, 
whereas the dashed line is the modulus of the interference term. 
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FIG. 5: The asymmetry term normalized to the symmetric terms as a function of 0. In this 
calculation the BFKL UGDF was used and —2 < yu < 2. The solid curve corresponds to the 
inclusive case while the dash-dotted curve is for the extra cut |p_l| > 50 GeV. 
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